Direct detection of WIMPs * 

David G. Cerdeno"^, Anne M. Green^ 

" Departamento de Fisica Teorica C-XI, and Institute de Fisica Teorica UAM-CSIC, 
Universidad Autonoma de Madrid, Cantoblanco, E-28049 Madrid, Spain 
School of Physics and Astronomy, University of Nottingham 
University Park, Nottingham, NG7 2RD, UK 

February 9, 2010 
Abstract 

A generic weakly interacting massive particle (WIMP) is one of the most attractive candidates 
to account for the cold dark matter in our Universe, since it would be thermally produced with 
the correct abundance to account for the observed dark matter density. WIMPs can be searched 
for directly through their elastic scattering with a target material, and a variety of experiments are 
currently operating or planned with this aim. In these notes we overview the theoretical calculation 
of the direct detection rate of WIMPs as well as the different detection signals. We discuss the 
various ingredients (from particle physics and astrophysics) that enter the calculation and review the 
theoretical predictions for the direct detection of WIMPs in particle physics models. 

1 Introduction 

If the Milky Way's DM halo is composed of WIMPs, then the WIMP flux on the Earth is of 
order 10^(100 GeV/m^) cm s . This flux is sufficiently large that, even though the WIMPs 
are weakly interacting, a small but potentially measurable fraction will elastically scatter off 
nuclei. Direct detection experiments aim to detect WIMPs via the nuclear recoils, caused 
by WIMP elastic scattering, in dedicated low background detectors [l[. More specifically 
they aim to measure the rate, R, and energies, Eft, of the nuclear recoils (and in directional 
experiments the directions as well). 

In this chapter we overview the theoretical calculation of the direct detection event rate 
and the potential direct detection signals. Sec. [2] outlines the calculation of the event rate, in- 
cluding the spin independent and dependent contributions and the hadronic matrix elements. 
Sec. [3] discusses the astrophysical input into the event rate calculation, including the local 
WIMP velocity distribution and density. In Sec.[l]we describe the direction detection signals, 
specifically the energy, time and direction dependence of the event rate. Finally in Sec. [5] we 
discuss the predicted ranges for the WIMP mass and cross-sections in various particle physics 
models. 



*This contribution appeared as chapter 17, pp. 347-369, of "Particle Dark Matter: Observations, Models 
and Searches" edited by Gianfranco Bertone, Copyright 2010 Cambridge University Press. Hardback ISBN 
9780521763684, |http://camb ridge.org/us/c atalogue/cata logue.asp?isbn=9780521 763684^ 
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2 Event rate 



The differential event rate, usually expressed in terms of counts/kg/day/keV (a quantity 
referred to as a differential rate unit or dru) for a WIMP with mass niy. and a nucleus with 
mass m]\f is given by 

vf{v) [v,En)dv, (1) 



dER rriN 7^„,„ dER 

where po is the local WIMP density, "'^^^ (f , Er) is the differential cross-section for the 
WIMP-nucleus elastic scattering and f{v) is the WIMP speed distribution in the detector 
frame normalized to unity. 

Since the WIMP-nucleon relative speed is of order 100km~^s~^ the elastic scattering 
occurs in the extreme non-relativistic limit, and the recoil energy of the nucleon is easily 
calculated in terms of the scattering angle in the center of mass frame, 9* 

^ //>^(i-cosr) ^ 

mN 

where pN = fn^rriN / {m-^ + rriN) is the WIMP-nucleus reduced mass. 

The lower limit of the integration over WIMP speeds is given by the minimum WIMP 



speed which can cause a recoil of energy Er: Vmin = y {n^N Er) / {2p'j^) . The upper limit is 
formally infinite, however the local escape speed Vesc (see Sec. 13. 2p . is the maximum speed in 
the Galactic rest frame for WIMPs which are gravitationally bound to the Milky Way. 

The total event rate (per kilogram per day) is found by integrating the differential event 
rate over all the possible recoil energies: 

^= r^^^;^ r -n-)'-s^i-^En)dv, (3) 

JEt rriN Jv^in dbjR 

where Et is the threshold energy, the smallest recoil energy which the detector is capable of 
measuring. 

The WIMP-nucleus differential cross section encodes the particle physics inputs (and as- 
sociated uncertainties) including the WIMP interaction properties. It depends fundamentally 
on the WIMP-quark interaction strength, which is calculated from the microscopic description 
of the model, in terms of an effective Lagrangian describing the interaction of the particular 
WIMP candidate with quarks and gluons. The resulting cross section is then promoted to 
a WIMP-nucleon cross section. This entails the use of hadronic matrix elements, which de- 
scribe the nucleon content in quarks and gluons, and are subject to large uncertainties. In 
general, the WIMP-nucleus cross section can be separated into a spin-independent (scalar) 
and a spin-dependent contribution, 

dawN _ f dawN \ f dfWN\ 

dER ~ \ dER J SI V dER Jsd' ^ ^ 

Finally, the total WIMP-nucleus cross section is calculated by adding coherently the above 
spin and scalar components, using nuclear wave functions. The form factor, F(Er), encodes 
the dependence on the momentum transfer, q = ^J2m^ER, and accounts for the coherence 
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loss which leads to a suppression in the event rate for heavy WIMPs or nucleons. In general, 
we can express the differential cross section as 



dawN _ mN 
dEn " 



{a^'Fl,{En) + a'o''FsDiER)) , (5) 



where a^^'^^ are the spin-independent and -dependent cross sections at zero momentum 
transfer. 

The origin of the different contributions is best understood at the microscopic level, by 
analysing the Lagrangian which describes the WIMP interactions with quarks. The contribu- 
tions to the spin- independent cross section arise from scalar and vector couplings to quarks, 
whereas the spin-dependent part of the cross section originates from axial-vector couplings. 
These contributions are characteristic of the particular WIMP candidate (see, e.g., |2|]) and 
can be potentially useful for their discrimination in direct detection experiments. 

2.1 Spin-dependent contribution 

The contributions to the spin-dependent (SD) part of the WIMP-nucleus scattering cross 
section arise from couplings of the WIMP field to the quark axial current, 97^75 For 
example, if the WIMP is a (Dirac or Majorana) fermion, such as the lightest neutralino in 
supersymmetric models, the Lagrangian can contain the term 

£ D ag(X7''75X) (97^759) • (6) 

If the WIMP is a spin 1 field, such as in the case of LKP and LTP, the interaction term is 
slightly different, 

C D ay'^P^Bp d, B,){qri5q) • (7) 
In both cases, the nucleus, A^, matrix element reads 

{N\q^^^5q\N) = 2X^ {N\Jn\N) , (8) 

where the coefficients relate the quark spin matrix elements to the angular momentum of 
the nucleons. They can be parametrized as 



where J is the total angular momentum of the nucleus, the quantities Ag" are related to 
t of 



the matrix element of the axial- vector current in a nucleon, 1 97/^75'? I = 2s^"^Ag"\ and 
{Sp^n) = (-^|5'p,„|A^) is the expectation value of the spin content of the proton or neutron 
group in the nucleuqj. Adding the contributions from the different quarks, it is customary to 
define 

«P= E a„= E -^A^ (10) 

q=u,d,s ^^^^ q=uAs ^^^^ 



^These quantities can be determined from simple nuclear models. For example, the single-particle shell 
model assumes the nuclear spin is solely due to the spin of the single unpaired proton or neutron, and therefore 
vanishes for even nuclei. More accurate results can be obtained by using detailed nuclear calculations. 
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and 

A = j[ap{Sp)+an{Sn)]. (11) 

The resulting differential cross section can then be expressed (in the case of a fermionic 
WIMP) as 

(^) =^A^G|J(J + l)^. (12) 

V dEn J SD TTv^ 'S'(O) 

(using = 2mj<[dEji). The expression for a spin 1 WIMP can be found, e.g., in Ref. Q]- 
In the parametrization of the form factor it is common to use a decomposition into 
isoscalar, oq = + a„, and isovector, ai = — an, couplings 

S{q) = alSQQ{q) + aoai5'oi(g) + a\Sii{q) , (13) 

where the parameters Sij are determined experimentally. 

2.2 Spin-independent contribution 

Spin-independent (SI) contributions to the total cross section may arise from scalar-scalar 
and vector-vector couplings in the Lagrangian: 

D Oi^XXm + agXlfiXQI^Q ■ (14) 

The presence of these couplings depends on the particle physics model underlying the WIMP 
candidate. In general one can write 

[dEnJsi~ 2/.>2 ' i^^i 

where the nuclear form factor for coherent interactions F'^{Er) can be qualitatively under- 
stood as a Fourier transform of the nucleon density and is usually parametrized in terms of 
the momentum transfer as [1; 0] 

where ji is a spherical Bessel function, s ~ 1 fm is a measure of the nuclear skin thickness, 
and Ri = \/ B? — 5s^ with R ~ 1.2^4^/^ fm. The form factor is normalized to unity at zero 
momentum transfer, -F(O) = 1. 

The contribution from the scalar coupling leads to the following expression for the WIMP- 
nucleon cross section, 

ao = ^[Z/f + (^-Z)rf , (17) 
vr 

with 

— = E — + E — ' (18) 

mj, m„ ^ 27 ^ rUn 

^ q=u,d,s ^ q=c,b,t ^ 

where the quantities fj^^ represent the contributions of the light quarks to the mass of the 
proton, and are defined as rupfj^^ = {p\mgqq\p). Similarly the second term is due to the 



4 



interaction of the WIMP and the gluon scalar density in the nucleon, with fj,Q = 1 — 
J2q=u,d,s frq- They are determined experimentally, 

/P^ = 0.020 ± 0.004, /^^ = 0.026 ± 0.005, /^^ = 0.118 ± 0.062, (19) 

with /ji^ = fj.^, fj,^ = fj.^, and /^^ = fj,^. The uncertainties in these quantities, among 
which the most important is that on f^s, mainly stem from the determination of the vr-nucleon 
sigma term. 

The vector coupling (which is present, for example, in the case of a Dirac fermion but 
vanishes for Majorana particles) gives rise to an extra contribution. Interestingly, the sea 
quarks and gluons do not contribute to the vector current. Only valence quarks contribute, 
leading to the following expression 

^° = -64^ ' ^2°) 

with 

BN = a^(.A + Z) + a^{2A-Z). (21) 

Thus, for a general WIMP with both scalar and vector interactions, the spin-independent 
contribution to the scattering cross section would read 



daw n\ 2 TTiTv 
dER J SI TTw 
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[zr + iA-z)rf + ^ 



F\Er) . (22) 



In most cases the WIMP coupling to neutrons and protons is very similar, w and 
therefore the scalar contribution can be approximated by 

(^) = ""''^T^^ HE.). (23) 

V dEn J SI TTv^ 

The spin-independent contribution basically scales as the square of the number of nucleons 
{A"^), whereas the spin-dependent one is proportional to a function of the nuclear angular 
momentum, (J + 1)/J. Although in general both have to be taken into account, the scalar 
component dominates for heavy targets {A > 20), which is the case for most experiments 
(usually based on targets with heavy nuclei such as Silicon, Germanium, Iodine or Xenon). 
Nevertheless, dedicated experiments exist that are also sensitive to the SD WIMP coupling 
through the choice of targets with a large nuclear angular momentum. 

As we have seen, the WIMP direct detection rate depends on both astrophysical input 
(the local DM density and velocity distribution, in the lab frame) and particle physics input 
(nuclear form factors and interaction cross-sections, which depend on the theoretical frame- 
work in which the WIMP candidate arises). We will discuss these inputs in more detail in 
Sees. [3] and [5] respectively. 



2.3 Hadronic Matrix Elements 

The effect of uncertainties in the hadronic matrix elements has been studied in detail for the 
specific case of neutralino dark matter ^5rt9i]. Concerning the SI cross section, the quantities 
fj^^ in Eq. (jl9p can beparametrized in terms of the vr nucleon sigma term, S^^r, (see in this 
respect, e.g., Refs.0; [^) which, in terms of the u and d quark masses reads 

^ttN = ^ (m„ + ma) {N\uu + dd\N) , (24) 
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and is related to the strange quark scalar density in the nucleon. The largest source of 
uncertainty in fj,^ stems from the determination of this quantity, for which the current 
data implies S^at = (64 ± 8) MeV [l^, which translates into a variation of a factor 4 in 
Jts- Notice that in general the WIMP interaction with strange quarks would be the leading 
contribution to the SI cross section, due to its larger Yukawa coupling. In this case, o"o is 
roughly proportional to fxs ^ and the above uncertainty in the strange quark content leads 
to a variation of more than one order of magnitude in the resulting SI cross section 0; @| ) • 

Similarly, for the SD cross section the uncertainties in the strange spin contribution 
are the dominant contribution to the error in gq. However, in the case of the neutralino, 
this can imply a correction of as much as a factor 2 in the resulting cross section 0], being 
therefore much smaller than the above uncertainty for the SI cross section. It should be 
emphasized, however, that uncertainties in the determination of the spin form factors S{q) 
would also affect the theoretical predictions for the dark matter detection rate. 



3 Astrophysics input 

3.1 Local DM density 

The differential event rate is dire ctly proportional to the local WIMP density, po = p{f = 
Rq) where Rq = (8.0 it 0.5) kpc 11[ is the solar radius. Any observational uncertainty in 
Pq therefore translates directly into in an uncertainty in the event rate and the inferred 
constraints on, or measurements of, the scattering cross-sections. 

Exclusion limits are traditionally calculated assuming a canonical local WIMP density, 
po = 0.3GeVcm~'^. The local WIMP density is calculated by applying observational con- 
straints (including measurements of the rotation curve) to models of the Milky Way and 
the values obtained can vary by a factor or order 2 depending on the models used [12Ml5l|. 
A recent study [16], using spherical halo models with a cusp {p oc r(r)~" as r — )• 0) finds 
po = (0.30 ± 0.05) GeVcm-3. 

3.2 Speed distribution 

The standard halo model, conventionally used in calculations of exclusion limits and signals, 
has an isotropic, Gaussian velocity distribution (often referred to as Maxwellian) 

= '''' 

The speed dispersion is related to the local circular speed by o" = and Vc = (220 ± 

20)kms-^ [17j] (see Sec JS.Sp so that a ~ 270 km s ^. This velocity distribution corresponds 
to an isotropic singular isothermal sphere with density profile p{r) oc r~^. The isothermal 
sphere is simple, and not unreasonable as a first approximation, however it is unlikely to be an 
accurate model of the actual density and velocity distribution of the Milky Way. Observations 
and numerical simulations (see chapter 1) indicate that dark matter halos do not have a 1/r^ 
density profile and are (to some extent at least) triaxial and anisotropic. 

If the velocity distribution is isotropic there is a one to one relation between /(v) and the 
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the spherically symmetric density profile given by Eddington's formula [l8^, see Refs. 
In general the steady state phase-space distribution of a collection of collisionless particles 
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is given by the collisionless Boltzmann equation and the velocity dispersions of the system 
are calculated via the Jean's equations (e.g. Ref. [2l[). Solving the Jean's equations requires 
assumptions to be made, and therefore even for a specific density distribution the solution is 
not unique. Several specific models have been used in the context of WIMP direct detection 



signals. The logarithmic ellipsoidal model 2j] is the simplest triaxial generalisation of the 



isothermal sphere and has a velocity distribution which is a multi-variate Gaussian. Osipkov- 
Merritt models 2^; 24 ] are spherically symmetric with radially dependent anisotropic velocity 



distributions. Fitting functions for the speed distributions in these models are available, for 
a selection of density profiles, in Ref. [25i]. Velocity distributions have also been extracted 



from cosmological simulations, with both multi-variate Gaussian |26l | and Tsallis j27l | distri- 
butions [2^ being advocated as fitting functions. While it is not known whether any of these 
models provide a good approximation to the real local velocity distribution function, the 
models are none the less useful for assessing the uncertainties in the direct detection signals. 

Particles with speed, in the Galactic rest frame, greater than the local escape speed, 
Vesc = V^I^C^oTl where $(r) is the potential, are not gravitationally bound. Many of the 
models used, in particular the standard halo model, formally extend to infinite radii and 
therefore their speed distribution has to be truncated at Vesc 'by hand' (see e.g. Ref. (29|). 
The standard value for the escape speed is Vesc = 650kms~^. A recent analysis, using high 
velocity stars from the RAVE survey, finds 498km s^^ < Vesc < 608kms~^ with a median 
likelihood of 544 km s"^ [30]. 

In Sec. m we discuss the impact of uncertainty in the speed distribution on the direct 
detection signals. 

3.3 Earth's motion 

The WIMP speed distribution in the detector rest frame is calculated by carrying out, a time 
dependent, Galilean transformation: v — t- v = v + Ve(t). The Earth's motion relative to the 
Galactic rest frame, Ve(t), is made up of three components: the motion of the Local Standard 
of Rest (LSR), the Sun's peculiar motion with respect to the LSR, Vq, and the Earth's orbit 
about the Sun, v°''^ . 

If the Milky Way is axisymmetric then the motion of the LSR is given by the local circular 
velocity {0,Vc,0), where Vc = 220kms~^ is the standard value. Kerr and Lynden-Bell found, 
by combining a large number of independent measurements, Vc = (222 it 20)kms~^ [l7 
more recent determination, using the proper motions of Cepheids measured by Hipparcos 
is broadly consistent: Vc = (218 ± 7) kms~"'^(-Ro/8kpc). 

The Sun's peculiar motion, determined using the parallaxes and proper motions of stars 
in the solar neighbourhood from the Hipparcos catalogue, is Vq = (10.0 it 0.4, 5.2 it 0.6, 7.2 it 
0.4) kms^"*^ [32] in Galactic co-ordinates (where x points towards the Galactic center, y is the 
direction of Galactic rotation and z towards the North Galactic Pole). 

A relatively simple, and reasonably accurate, expression for the Earth's motion about 
the Sun can be found by ignoring the ellipticity of the Earth's orbit and the non-uniform 
motion of the Sun in right ascension [s^: v"""^ = Ve[eisinX{t) — e2CosA(t)] where Ve = 
29.8 km s~^ is the orbital speed of the Earth, X{t) = 27r(t — 0.218) is the Sun's ecliptic 
longitude (with t in years) and ex(2) are unit vectors in the direction of the Sun at the Spring 
equinox (Summer solstice). In Galactic co-ordinates ei = (—0.0670,0.4927,-0.8676) and 



31|, 



62 = (—0.9931, —0.1170,0.01032). More accurate expressions can be found in Ref. 341 ]. 
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The main characteristics of the WIMP signals can be found using only the motion of 
the LSR, and for the time dependence the component of the Earth's orbital velocity in that 
direction. However accurate calculations, for instance for comparison with data, require all 
the components described above to be taken into account. 



3.4 Ultra-fine structure 

Most of the WIMP velocity distributions discussed in Sec. 13.21 are derived by solving the 
coUisionless Boltzmann equation, which assumes that the phase space distribution has reached 
a steady state. However this may not be a good assumption for the Milky Way; structure 
formation in CDM cosmologies occurs hierarchically and the relevant dynamical timescales 
for the Milky Way are not many orders of magnitude smaller than the age of the Universe. 

Both astronomical observations and numerical simulations (due to their finite resolution) 
typically probe the dark matter distribution on ~ kpc scales. Direct detection experiments 
probe the DM distribution on sub milli-pc scales (the Earth's speed with respect to the 



Galactic rest frame is ^ 0.2mpcyr~^). It has been argued (see e.g. Refs. |35l437l]) that on 
these scales the DM may not have yet reached a steady state and could have a non-smooth 
phase-space distribution. On the other hand it has been argued that the rapid decrease in 
density of streams evolving in a realistic, ellipsoidal. Galactic potential means that there will 
be a large number of over lapp ing streams in the solar neighbourhood producing an effectively 



smooth DM distribution 381]. 



If the local DM distribution consists of a small number of streams, rather than a smooth 
distribution, then there would be significant changes in the signals which we will discuss in 
Sec. m This is currently an open issue; directly calculating the DM distribution on the scales 
probed by direct detection experiments is a difficult and unresolved problem. 

It has been suggested that a tidal stream from the Sagittarius (Sgr) dwarf galaxy, which 
is in the process of being disrupted, passes through the solar neighbourhood with the asso- 
ciated DM potentially producing distinctive signals in direct detection experiments [iol: l4o|. 
Subsequent numerical simulations of the disruption of Sgr along with observational searches 
for local streams of stars suggest that the Sgr stream does not in fact pass through the solar 



neighbourhood (see e.g . Ref. 41] )• None the less the calculations of the resulting WIMP 
signals in Refs. 3^; 4C| provide a useful illustration of the qualitative effects of streams. 



4 Signals 

We have already seen in Sec. [2] that the recoil rate is energy dependent due to the kinematics 
of elastic scattering, combined with the WIMP speed distribution. Due to the motion of the 
Earth with respect to the Galactic rest frame the recoil rate is also both time and direction 
dependent. In this section we examine the energy, time and direction dependence of the recoil 
rate and the resulting WIMP signals. In each case we first focus on the signal expected for 
the standard halo model, with a Maxwellian velocity distribution, before discussing the effect 
on the signal of changes in the WIMP velocity distribution. 
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4.1 Energy dependence 

The shape of the differential event rate depends on the WIMP and target masses, the WIMP 
velocity distribution and the form factor. For the standard halo model the expression for the 



differential event rate, eq. [H can be rewritten approximately (c.f. Ref. 4J]) as 

dR f dR\ ^o,^ , ( Er\ , , 

where {dRjdE^Q is the event rate in the -E — t- keV limit. The characteristic energy scale 
is given by E^ = (ci2/i^t;^)/mjv where ci is a parameter of order unity which depends on 
the target nuclei. If the WIMP is much lighter than the target nuclei, rriy. <C m^, then 
Ec oc /mj\j- while if the WIMP is much heavier than the target nuclei E^ oc mjv . The total 
recoil rate is directly proportional to the WIMP number density, which varies as ^/m^. 

In fig. [T] we plot the differential event rate for Ge and Xe targets and a range of WIMP 
masses. As expected, for a fixed target the differential event rate decreases more rapidly with 
increasing recoil energy for light WIMPs. For a fixed WIMP mass the decline of the differen- 
tial event rate is steepest for heavy target nuclei. The dependence of the energy spectrum on 
the WIMP mass allows the WIMP mass to be estimated from the energies of detected events 
(e.g. Ref. [i^). Furthermore the consistency of energy spectra measured by experiments 
using different target nuclei would confirm that the events were due to WIMP scattering 
(rather than, for instance, neutron backgrounds) In particular, for spin independent 

interactions, the total event rate scales as A"^. The is sometimes referred to as the 'materials 
signal'. 

The WIMP and target mass dependence of the differential event rate also have some 
general consequences for experiments. The dependence of the total event rate on rriy. means 
that, for fixed cross-section, a larger target mass will be required to detect heavy WIMPs 
than lighter WIMPs. For very light WIMPs the rapid decrease of the energy spectrum with 
increasing recoil energy means that the event rate above the detector threshold energy, E'j' , 
may be small. If the WIMP is light, < 0(10 GeV), a detector with a low, < ©(keV), 
threshold energy will be required. 

The most significant astrophysical uncertainties in the differential event rate come from 
the uncertainties in the local WIMP density and circular velocity. As discussed in Sec. 13. li the 
uncertainty in the local DM density translates directly into an uncertainty in constraints on 
(or in the future measurements of) the scattering cross-section. The time aiiera^ed differential 
event rate is found by integrating the WIMP velocity distribution, therefore it is only weakly 
sensitive to changes in the shape of the WIMP velocity distribution. For the smooth halo 
models discussed in Sec. 13.21 the time averaged differential event rates are fairly similar to 
that produced by the standard halo model [3; . Consequently exclusion limits vary only 



weakly 45|; 140] and there would be a small (of order a few per-cent) systematic uncertainty 



in the WIMP mass deduced from a measured energy spectrum 471]. With multiple detectors 
it would in principle be possible to measure the WIMP mass without any assumptions about 
the WIMP velocity distribution [48]. 

In the extreme case of the WIMP distribution being composed of a small number of 
streams the differential event rate would consists of a series of (sloping due to the form 
factor) steps. The positions of the steps would depend on the stream velocities and the 
target and WIMP masses, while the relative heights of the steps would depend on the stream 
densities. 
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Figure 1: The dependence of the spin independent differential event rate on the WIMP mass 
and target. The sohd and dashed hnes are for Ge and Xe respectively and WIMP masses of 
(from top to bottom at Er = OkeV) 50, 100 and 200 keV. The scattering cross-section on 
the proton is taken to be dp ^ = 10~* pb. 



4.2 Time dependence 

The Earth's orbit about the Sun leads to a time dependence, specifically an annual modula- 



tion, in the differential event rate [29l : |49| |. The Earth's speed with respect to the Galactic 
rest frame is largest in Summer when the component of the Earth's orbital velocity in the 
direction of solar motion is largest. Therefore the number of WIMPs with high (low) speeds 
in the detector rest frame is largest (smallest) in Summer. Consequently the differential event 
rate has an annual modulation, with a peak in Winter for small recoil energies and in Summer 
for larger recoil energies [50|]- The energy at which the annual modulation changes phase is 
often referred to as the 'crossing energy'. 

Since the Earth's orbital speed is significantly smaller than the Sun's circular speed the 
amplitude of the modulation is small and, to a first approximation, the differential event rate 
can, for the standard halo model, be written approximately as a Taylor series: 

[1 + A{Er) cos a{t)] , (27) 



dR ^ / dR\ 



where a{t) = 27r(t — tQ)/T, T = 1 year and Iq ~ 150 days. In fig. [2] we plot the energy 

1/2 

dependence of the amplitude in terms of f min (recall that Vmin oa E^ with the constant of 
proportionality depending on the WIMP and target nuclei masses). The amplitude of the 
modulation is of order 1-10 %. 

The Earth's rotation provides another potential time dependence in the form of a diur- 



nal modulation as the Earth acts as a shield in front of the detector 5l|; l52], however the 
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Figure 2: Dependence of the amplitude of the annual modulation, A(?7min)i on Vr, 



amplitude of this effect is expected to be small, < 1% j52[ |. 

There has been a substantial amount of work on the annual modulation for the non- 
Maxwellian velocity distributions described in Sec. 13.21 [l^ : [20I : 53-59]. In contrast to the 



time-averaged differential event rate, both the phase and amplitude of the annual modulation 
can vary substantially. Consequently the regions of WIMP mass-cross-section parameter 
space consistent with an observed signal can change significantly 5J; |60(]. Note that if the 
components of the Earth's orbital velocity perpendicular to the direction of Solar motion are 
neglected, then the phase change will be missed 53]. For a WIMP stream the position and 
height of the step in the energy spectrum would vary annually (e.g. Ref. [giI]). 



4.3 Direction dependence 

The detector motion with respect to the Galactic rest frame also produces a directional signal. 
The WIMP flux in the lab frame is sharply peaked in the direction of motion of the Earth, 
and hence the recoil spectrum is also peaked in this direction (albeit less sharply due to the 
elastic scattering). 

The directional recoil rate is most compactly written as [s^j 



dR 



dEndn 



(28) 



where dQ = d(j)dcos9, q is the recoil direction and f{vmm,^) is the 3-dimensional Radon 
transform of the WIMP velocity distribution /(v) 



f{Vn 



J(v.q 



(29) 
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Geometrically the Radon transform is the integral of the function / (v) on a plane orthogonal 
to the direction q at a distance Vmm from the origin. See Ref. [63.] for an alternative, but 
equivalent, expression. 

For the standard halo model the direction dependence is approximately given by [6j] 



dR 



dER d cos 7 



cx exp 



{vq cos 7 



(30) 



where 7 is the angle between the recoil and the mean direction of solar motion. The dis- 
tribution of recoil directions peaks in the mean direction of motion of the Sun (towards the 
constellation CYGNUS ^]) with the event rate in the forward direction being roughly an 
order of magnitude larger than that in the backward direction [6j], since the Sun's speed is 
comparable to the mean WIMP speed. The directional recoil rate of = 100 GeV WIMPs 
on S is shown in fig. El 

With an ideal directional detector, with 3-d read-out and capable of measuring the senses 
of the recoils {i.e. distinguishing between the head and tail of each recoil), it would be 
possible to distinguish a WIMP signal from isotropic backgrounds with only of order 10 
events 6^; 6^; 63]. This number increases significantly (by roughly an order of magnitude) if 
either the senses can not be measured or the read out is only 2-d Another potential 

directional signal is the rotation of the mean recoil direction in the lab over a sidereal day due 
to the motion of the Earth [63] • See chapter 22 for the principles and practice of directional 
detection experiments. 

For plausible smooth halo models changes in the WIMP velocity distribution affect the 
detailed angular recoil rate. However the rear-front asymmetry is robust and the number of 
events required by an ideal detector only varies by of order 10% [62 



66l: 16711 . For non-ideal 



detectors the variation in the number of events required can be larger 67|; 1681]. With a large 
number of events (of order thousands) it would be possible to probe the detailed WIMP 



velocity distribution 67l; l68|; l7l|| . A stream of WIMPs produces a recoil spectrum which is 



peaked in the opposite direction (e.g. Ref. |72l|). 



5 Particle Physics input 

Let us finally address the Particle Physics input to the determination of the WIMP detec- 
tion rate, which enters through the theoretical predictions to the WIMP-nucleus scattering 
cross section. These are sensitive to the WIMP nature. We will here briefiy summarise the 
results for WIMP candidates that arise in various well-motivated theories of physics beyond 
the Standard Model at the TeV scale (Supersymmetric theories, Universal Extra Dimension 
scenarios and Little Higgs models), as well as phenomenologically motivated scenarios. 



5.1 Supersymmetric WIMPs 

The canonical and best studied supersymmetric WIMP is the lightest neutralino, Xi- Its 
detection properties are very dependent on its composition. More specifically, within the 
MSSM framework, in the expressions for the scattering amplitudes @; 173-76], the SI part of 
the neutralino-nucleon cross section receives contributions from Higgs exchange in a t-channel 
and squark exchange in an s-channel. The latter also contributes to the SD part of the cross 
section, together with a Z boson exchange in a t-channel. 



12 



+90° 




0.002 0.004 0.006 0.008 0.01 

Recoil Rate(E^>20ke V)/kg" day 'sr- ' 



Figure 3: Hammer- Aitoff projection of the directional recoil rate for S, above = 20keV, 
due to a standard halo of WIMPs with = 100 GeV. This figure was generated using the 
HADES directional dark matter simulation code written by Ben Morgan. 



The dependence of the cross sections, and detection prospects, on the neutralino compo- 
sition are well known. For example, a large Higgsino component induces an enhancement of 
both the Higgs and Z boson exchange diagrams, thereby leading to an increase in both the 
SD and SI cross sections. On the other hand, the presence of light squarks (if they become 
almost degenerate with the neutralino) can lead to an enhancement of (mainly) the SD cross 
section. 

Analyses of general super symmetric scenarios with parameters defined at low energy reveal 
that the neutralino SI cross section can be as large as 10~^ pb for a wide range of neutralino 
masses up to 1 TeV [S; 77]. Interestingly, when gaugino masses not fulfilling the GUT relation 
are allowed, very light neutralinos with masses m-^o > 7 GeV, and potentially large cross- 
sections, can be obtained 



78l : |82| |. It has been argued that these neutralinos could account 
for the DAMA/LIBRA annual modulation signal without contradicting the null results from 
CDMS and XENONIO [t^, however this interpretation is now more constrained 80]. All 
these features are illustrated in in FigUl 

Analyses have also been done from the point of view of Supergravity theories, where 
the SUSY parameters are defined at the GUT scale. In the simplest case, the CMSSM, 
the cross-sections are generally small, since Xi is mostly bino. The largest cross-sections, 
a^^ « 0(10"^ pb), are obtained in the focus point region, where the neutralino becomes a 



mixed bino-Higgsino state fss!; [sj]. Interestingly, this region seems to be favoured by recent 
Bayesian analyses of the CMSSM parameter space sB-Sa]. Moreover, the predicted SD cross 
section is also sizable in the focus point region, approximately reaching 10~^ — 10~^ pb. In 
more general Supergravity scenarios the predicted cross-sections can be signific antly larger 
through the inclusion of non-universal values for either the scalar masses [Sgl-Eoi] (non- 
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Figure 4: Left) Theoretical predictions for neutralino-nucleon SI cross section as a function of 
the neutrahno mass obtained by combining the scans of the MSSM parameters from Refs 



77h79|]. The theoretical predictions for the SI cross section of very light neutralinos in the 



NMSSM from Ref. {80i| are shown by means of empty grey circles. Present and projected 
experimental sensitivities are displayed using solid and dashed lines respectively. Right) 
Theoretical predictions for neutralino-nucleon SD cross section as a function of the neutralino 
mass, using the data from the supergravity scan of [8l|]. 



universalities in the Hi ggs mass parameters being the most effective), the gaugino masses 

Tooli Eosl - Eogl i. or both [ii3; [nil. 



The detection prospects of the lightest neutralino in extended supersymmetric models 
may be significantly different, mostly due to the changes in the Higgs sector and the presence 
of new neutralino states. The se constructions are generally referred to as singlet extensions 



of the MSSM (see, e.g., |112l |) For example, this is the case of the NMSSM, in which the 



presence of very light Higgses (consistent with LEP constraints if t hey have a large singlet 
composition) can lead to a sizable increase of the SI cross section jll3 - 115 |. Mo reover, in 
the NMSSM very light neutralinos (with masses below 10 GeV) are viable 1161 ] and can 
have very distinctive predictions for their direct detection, including, for example, smaller SI 
cross section than in the MSSM [scJ]. The theoretical predictions for the SI cross section of 
neutralinos with m o < 30 GeV are plotted in FigHl In general, the singlet component of the 
neutralino does not couple to the Z boson or to squarks and thus in these constructions the 
theoretical predictions for the SD cross section remain the same as in the MSSM. 

Finally, there is another viable supersymmetric WIMP candidate for dark matter, the 
lightest sneutrino. The left-handed sneutrino in the MSSM is excluded given its sizable 
coupling to the Z boson. They therefore either annihilate too rapidly, resulting in a very small 
relic abundance, or have larg e sca ttering cross sections which have already been excluded by 
direct detection experiments 117l |. Several models have been proposed to revive sneutrino DM 
by reducing its coupling with the Z-boson . This can be achieved by introducing a mixture of 
left- and right-handed sneutrinos 118l4l2ll | , or by considering a purely right-handed sneutrino 
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in m odels with an extended gauge sector 122l | or Higgs sector 123 : 124 1 such as the NMSSM 



125l |. In the first class of models, the elastic scattering of sneutrinos with quarks would take 
place through the i-exchange of Z bosons, whereas in the second class it would mostly be due 
to the exchange of Higgs bosons. The resulting SI cross section in these cas es can be within 
the reach of future detectors for a wide range of sneutrino masses |124l : 125l |. Being a scalar 
particle, the SD cross section vanishes for the sneutrino. 

5.2 Kaluza Klein dark matter in UED 

Mode ls of Universal Extra Dimensions, in which all fields are allowed to propagate in the 
bulk [liet . also provide well-motivated candidates for WIMP dark matter in the form of 



the Lightest Kaluza-Klein Particle (LKP), which is likely to be associated with the first KK 
excitation of the hypercharge gauge boson, -BJ^^ [l27l:[l28l|. 



The elastic scattering of b'^^ with quarks takes place through the exchange of KK quarks 
along t and s channels, which contribute to both the SD and SI cross sec tion, and a Higgs 
exchange along a t-channel which only gives a SI contribution 12^131]. The theoretical 



predictions for the elastic scattering cross sections of are very dependent on the mass 
splitting between it and the KK quark, Ag. In particular, both the SI and SD contributions 
increase when Ag becomes small, as a consequence of the enhancement of the contribution 
from KK quarks. This is relevant, since in the UED scenario one expects a quasi-degenerate 
spec trum, in which the splittings between the masses are only induced by radiative corrections 



1271] . The SI cross section can also be larger in the presence of light Higgses and for small LKP 
masses. However the Higgs mass is generally larger than in the MSSM and its contribution 
is suppressed with respect to that of KK quarks. 

It has been shown that the theoretical predictions for the SI cross section of B^^ can be 
as large as a^^ ~ 10~^ pb for masses ranging fr om 500 G eV to 1 TeV, when Ag ~ 0.01, for 
which the correct relic density can be obtained HH; flii ]. Under the same conditions, the 



predicted SD cross section are smaller and, for masses up to 1 TeV ton-scale detectors would 



be required to detect them. These predictions are illustrated in Fig.[5l from [134]. 

Other LKP candidates are possible within the UED model if nonvanishing boundary 
terms are allowed. More specifically, one may c onsid er the first excited states associated 
with either the Z boson, Z(^), or the Higgs, H^^^ [isj- 

The detection properties of the Z^"^^ 
are very similar to those of the B^^ [134*] (although for the Z^^^ the neutron and proton 

spin-dependent cross section are exactly the same, contrary to the case of the Sj^^). 

Finally, in models with two universal extra dimensions the LKP generally corresponds 
to the KK e xcitation associated with the hypercharge gauge boson, which is called spinless 
photon [135]. Being a scalar, this particle has no SD cross section. Its SI cross section can 
be similar to that of the b'^^ . 

5.3 Little Higgs dark matter 

In these constructions a discrete symmetry, called T-parity, is introduced in order to allevi- 
ate the stringent experimental constraints on low-energy observables. A phenomenological 
consequence of T-parity is that the Lightest T-odd Particle (LTP) becomes absolut ely stable . 
Interestingly, the LTP is usually the partner of the hypercharge gauge boson Bh [l36l - [l38j . 



15 



10" 



10" 



10" 



g 10-^ 



10" 



10 



-10 



10 



-11 




. . ton-scale detector 



_i_ 



_i_ 





103 




10^ 




loi 




10° 


■3" 


10-1 


N y 


10~2 


a 
o 

o 


10-3 


o. 

.3 


10-4 


p. 

m 

b 


10-5 




10-6 




10-"? 




10-8 




10-9 



200 400 



600 800 
JCP (GeV) 



1000 1200 



XENONIO 



COUPP 




KIMS 

ton-scale detector 
7i LKP 



(b) 



Zi LKP 



200 400 600 800 
niLKp (GeV) 



1000 1200 



Figure 5: Left) Theoretical predictions for the SI LKP-nucleon scattering cross section as a 
function of the LKP mass for the LKP candidates discussed i n th e text. Right) Predictions 
for the SD LKP-proton cross section. Figures extracted from 134l |. 



The LTP-nucleon scattering cross section receives SI contributions via H iggs and heavy quark 
exchange, the latter being the only contribution to the SD part 139l |. Thus, the resulting 
expressions are very similar to the case of KK dark matter in UED. However, unlike the UED 
scenario, there is no reason for the heavy quarks to be degenerate in mass with Ah- This, 
together with the smallness of the heavy quark Yukawa couplings implies that their scatter- 
ing cross-sections are in fact very suppressed. The SI cross section, being dominated by the 
Higgs exchange t-channel, increases slightly when the Higgs mass is small but is nevertheless 
generally below IQ-^'^ pb. The theoretical predictions for the SD cross section are also very 
small. In summary, the direct detection of little Higgs dark matter is much more difficult 
than the SUSY and UED cases. 



5.4 Minimal Models and other approaches for dark matter 

Instead of considering DM candidates arising in existing theories beyond the SM, a bottom-up 
approach can be adopted in which minimal additions to the SM are considered, involving the 
inclusion of a WIMP field (usually a new singlet) and new symmetries that protect their decay 
(in some cases, also a new "mediator" sector that coupl es the WI MP to the SM). Examples 
in this direction include WIMPs with sing let mediation [l4ol - E4i |. models with an extended 
electroweak sect or |143l4l45| |. models with additional gauge groups, and the Secluded Dark 
Matter scenario |146l | in which WIMPs could escape direct detection. 

The theoretical predictions for the direct detection of WIMPs in this class of models 
are very dependent on the mediator sector, since it determines the couplings of WIMPs to 
ordinary matter (quarks). For example, scala r WIMPs interact with ordinary quarks through 
the exchange of Higgs bosons in a t channel 140l : 141]. The singlet coupling to the Higgs is 
constrained in order to reproduce the correct relic abundance, thus leaving only the WIMP 
and Higgs masses as free parameters. The resulting cross-sections increase as both masses 
decrease and can be as large as IQ-^ pb for very light WIMPs of order 10 GeV but are reduced 
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to be below 10~^ pb when the DM candidate has a mass above 100 GeV. 

Other scenarios can have more general couplings. In the Inert Doublet Model, where 
elastic scattering proceeds only through the exchange of a Higgs or a Z boson along a t 
channel, the resulting SI cross section is rather small. Only for light WIMPs, with masses 
below 100 GeV, is the predicted SI cross section large enough to be experimentally tested 
(from about 10~^'^pb to as much as 10~^pb) whereas the predictions for he avy WIMPs are 
well below the sensitivity of ton-size expe riments, usually of order 10~^^ pb [145]. 

In the Minimal DM approach of 143| | WIMP candidates with direct couplings to the Z 



boson are already excluded by direct DM searches. However, some fermionic candidates are 
still viable which interact with quarks through the exchange of W (and Higgs) bosons. These 
particles (with masses of several TeV) can have an SI cross section of order 10~^ pb. 

5.5 Inelastic cross section 

It is finally worth mentioning that the WIMP-nucleon cross section can also receive a con- 



tribution from inelastic scattering by creating either an excited ri uclear 147f | or electronic 



148f | or even by creating an excited WIMP state 149-151]. The last possibility is 



state 

particularly interesting if the mass difference, S, between the excited dark matter candidate, 
X2 and WIMP xi is of order of the WIMP kinetic energy (i.e., about 100 keV). In that case, 
the inelastic scattering off nuclei Xi^ ~^ X2^ can occur and the only kinematic change is in 
the minimal WIMP velocity that can trigger a specific recoil energy, which is increased by 
^Vmin = S / \/(^^tt^nEr) . This clearly favours detection in heavy targets such as Iodine (since 
^Vmin is Smaller) and might provide a possible explanation for the DAMA/LIBRA signal 
compatible with the null restuls in other experiments [1511] . 

5.6 Discrimination of Dark Matter candidates 

As illustrated by figs. U] and O current experiments are already probing the masses and 
cross-sections predicted for various WIMP candidates. Furthermore future experiments will 
be sensitive to a substantial fraction of the parameter space. If any of these experiments 
succeed in detecting dark matter particles, the next objective will be to identify its particle 
nature. 

In this sense, the simultaneous measurement of both the SI and SD dark matter couplings, 
through experiments which are sensitive to both signals, can provide very valuable informa- 
tion 8l[ . This is illustrated in Fig. [6l where the ratio of SI to SD cross section is plotted for 
the neutralino and LKP cases versus the event rate for two complementary choices of target 
materials (that could be used in the COUPP experiment). As shown in the left panel, the 
measurement of an event rate in a single detector does reduce the number of allowed models, 
but does not generally place significant constraints on coupling parameters or on the nature 
of the dark matter detected. However, as shown in the right panel, a subsequent detection 
using a second complementary target does substantially reduce the allowed range of coupling 
parameters, and allows, in most cases, an effective discrimination between neutralino and 
LKP candidates. 

This analysis can be extended to other dark matter candidates. Moreover, in order to 
eliminate the large astrophysical and theoretical uncertainties which affect t he da rk matter 
rates, the ratio of WIMP-proton and WIMP-neutron amplitudes can be used and com- 



pared for different target materials. For example, the comparative study of the ratios of SD 
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Figure 6: Left) The detection -with a given target such as CF3I can only loosely constrain 
models for the dark matter (blue squares for neutralinos, green circles for the LKP) in the 
ap^/Up^ versus count-rate plane. Red (magenta) dots show the many models consistent with 
a measurement of ~ 10~^ (10~^) counts/kg day. Right) measurement of the event rate in a 
second target such as C4F10, with lower sensitivity to spin-independent couplings, effectively 
reduces the remaining number of allowed models-orange (aqua) dots-and generally allows 
discrimination between the neutralino and the LKP. Figures extracted from 81 1. 



neutron to proton amplitudes can provide a good discrimination of dark matter models by 
distinguishing candidates for which the SD cross section is dominated by Z boson exchange 
(such as the neutralino in some regions of the parameter space) from those where the domi- 
nant channel is squark or KK quark exchange (such as the LKP or LTP). A similar analysis 
for the SI neutron to proton ratio can be used to disentangle models with dominant Higgs 
or Z boson exchange, however in practise the different target materials are less sensitive to 
these differences. Finally, the mass determination techniques described in the Sec. 14.11 can 
provide complementary information that could lead to more effective discrimination between 
the various dark matter models. 
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